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Deuteron radial moments for renormalized chiral potentials 
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Abstract. We calculate deuteron positive and negative radial moments involving any bilinear function of 
the deuteron S and D wave functions for renormalized OPE and TPE chiral potentials. The role played 
by the strong singularities of the potentials at the origin and the short distance insensitivity of the results 
when the potentials are fully iterated is emphasized as compared to realistic potentials. 
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Chiral dynamics has played an important role in the 
theoretical description of low energy hadronic reactions 1 
and so far is the only known vestige of the underlying 
fundamental QCD theory of strong interactions in nu- 
clear physics. There is a number of low energy theorems 
based on chiral symmetry which provide a quantitative 
and model independent insight into low energy processes 
involving pions and nucleons, due to the clear scale sepa- 
ration between nuclear physics and QCD. For compound 
systems which at low energies disclose their composite na- 
ture the theoretical description necessarily becomes very 
involved and probably dependent on arbitrary assump- 
tions. On the contrary, for weakly bound systems such as 
the deuterium nucleus one expects important simplifica- 
tions leading to a more scheme independent and possi- 
bly systematic description of these systems. This possibil- 
ity motivated the introduction of Effective Field Theory 
(EFT) approaches for nuclear physics based on the chi- 
ral symmetry of QCD, and the derivation of low energy 
theorems, as, for example, pion-deuteron scattering [3] (for 
comprehensive reviews see e.g. Ref. 0.5 ,6 ). In many cases 
most of the information needed for reactions involving the 
deuteron can be encoded by simple deuteron matrix ele- 
ments. 

Guided by earlier work 0IH], we have proposed 
I10U1 111 12111. 'ij to renormalize the NN interaction in a non- 
perturbative way, highlighting model independent long dis- 
tance correlations among physical observables. In our ap- 
proach the long distance chiral NN One Pion Exchange 
(OPE) and Two Pion Exchange (TPE) potentials, com- 
puted within perturbation theory in Refs. [T4lll5lll6| . are 
iterated to all orders in the Schrodinger equation very 
much in the spirit of the original Weinberg approach 0. 
However, some subtleties are found |l()U12lll7j . which im- 
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pose strong constraints on the admisible short distance 
physics based on orthogonality, uniqueness and finiteness 
of the results. 

In the 3 Si — 3 £q channel, the relative proton-neutron 
state for negative energy is described by the coupled equa- 
tions 



h + U s {r) 
U sd {r) 



d 

dr' 2 



U sd (r) 
% + U d (r 



+ 




= -7 



1) 



Here 7 = V M B, with B the deuteron binding energy and 
M the nucleon mass, U (r) = MV(r) are the reduced po- 
tentials and u(r) and w(r) are S- and D-wave deuteron re- 
duced wave functions respectively. At long distances they 
satisfy, 
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where 77 is the asymptotic D/S ratio parameter and As 
is the asymptotic normalization factor, which is such that 
the deuteron wave functions are normalized to unity. For 
conventions and numerical values of parameters we use 
Ref. PTUT^] throughout. 

In this work we report on the radial moments 
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for — 3 < n < 2 which appear in many situations of in- 
terest, such as the calculation of the matter radius, the 
deuteron quadrupole moment and deuteron magnetic mo- 
ment for the positive powers, as well as ird and Kd elastic 
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Fig. 1. The OPE and TPE deuteron wave functions, u (left) and w (right), as a function of the distance (in fm) compared to 
the Nijmegen II wave functions 1221. The asymptotic normalization u — > e _7r has been adopted and the asymptotic D/S ratio 
is taken i) = 0.0256(4) in the TPE case (for OPE -q = 0.026333). We use the set IV of chiral couplings (see Ref. [T2]). 



scattering and neutral pion photoproduction, jd — > 7r°(i, 
in the case of the negative powers. 

An important issue is the finiteness of the negative ra- 
dial moments, a topic which has been recently discussed 
for OPE [HltlMTT?] and TPE [23]. The remarkable find- 
ing is that chiral potentials |14U15U16| . when fully iterated, 
have an increasing number of finite inverse radial moments 
due to the near the origin singularities of the potential. 
They smoothen the short distance behaviour of the wave 
functions and hence improve the convergence of the in- 
verse radial moments. This is in sharp contrast with per- 
turbative approaches [2U , for which the perturbative wave 
functions diverge at the origin |11II12| 1 , or conventional 
(regular) phenomenological potentials [22] where the S- 
and D-wave short distance behaviour of the wave func- 
tions, u ~ r and w ~ r 5 respectively, is enough to render 
(l/r) u and (l/r 2 ) u finite, but produce divergent higher 
inverse moments. We illustrate the situation below. 

At distances much shorter than the pion Compton 
wavelength, the OPE potential behaves as 



complete analysis see Ref. 52)- The elimination of the di- 
verging exponential fixes 770PE = 0.0263. From this short 
distance behaviour of the wave functions, one finds that 
the {l/r) u and (l/r 2 ) u moments are finite for the OPE 
potential, while (l/r 3 )„ and higher moments diverge, as 
it would happen for a regular potential. 

The short distance behaviour of the TPE (NNLO) has 
been exploited in Ref. J2]- The potential at short dis- 
tances behaves as [P4"lll5lllfi| 
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with R s = 0, R sd = 2V2R, R d = 4i? and R = 3g 2 A M/32nf 2 
(= 1.07764 fm). This behaviour of the potential is strong 
enough to overcome the centrifugal barrier at short dis- 
tances, thus modifying the usual short distance behaviour 
of the wave functions, which can schematically be written 
as 
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and Cj = Mci are the low energy chiral couplings appear- 
ing in ttN scattering. As in the OPE case, this potential 
is strong enough at short distances to modify the short 
distance behaviour of the wave function, which now reads 

u(r) ~ w(r) ~ C + (r/R + f 2 /+ (r/R + ) + 

C_ (r/R_f 2 f-(r/R-) (10) 
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where R+ and R_ are the eigenvalues of the matrix in 
Eq. ©, and f±(r/R±) represents a linear combination 
of sin {R\/2r 2 ) and cos {R\/2r 2 ) leaving ^tpe as a free 
parameter |12U2l)j . From this short distance behaviour, 

(4 V ^), 'cos(4 v /i!A) and cxp {-4*/2y/Rp) (for a the (V^m (V^ 2 )^ and (l/r 3 ) u radial moments are fi- 

" nite for the TPE potential, while higher moments diverge 

(although they would become finite for higher order po- 
tentials). 



where f(r/R) represents some linear combination of 



1 This divergency holds for perturbations both on boundary 
conditions or on distorded (fully iterated) OPE waves. 
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Table 1. Deuteron radial moments (in units of powers of fm). We consider the OPE and TPE potentials; in the case of the 
OPE potential we have taken g^NN = 13.08 (i.e. gA = 1.29, OPE) and gA = 1-26 (OPE*), while in the TPE case we show the 
results corresponding to the four set of chiral couplings considered along our previous works 12 20 . In the OPE case the error is 
estimated by varying the semiclassical matching radius |11I20) in the 0.1 — 0.2 fm range , while in the TPE case the error comes 
from the experimental uncertainty of the D/S ratio, r\ — 0.0256(4). TPE Sets I,II,II and IV refer to the chiral parameters, ci, 
C3 and C4 of Refs. |26) . |16) . |24) and |25| respectively. Nijm II and Reid93 are calculated from Ref. [22] or taken from Ref. |27| . 



The wave functions for OPE and TPE as compared to 
the Nijmll ones have been depicted in Fig. ^ The radial 
moments are tabulated in Tabled As we see, and despite 
the very different behaviour at short distances between 
the deuteron wave functions corresponding to renormal- 
ized chiral potentials and to phenomenological potentials, 
the convergent moments are fairly similar (the u+w com- 
bination works better) despite that Nijmll and Reid93 
contain no explicit TPE components. For inverse moments 
the trend improves clearly when going from OPE to TPE, 
which we interpret as a correct implementation of model 
independent long distance correlations generated by chiral 
symmetry and renormalization constraints. 
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